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Abstract 

We discuss the three independent asymmetries, A^\q 2 ), A^ m \q 2 ) and A^ e \q 2 ), 
that one can build from the amplitudes A±(q 2 ) and A\\(q 2 ). These quantities are 
expected to be accessible from the new B- physics experiments, they are sensitive to 
the presence of new physics, and they are not very sensitive to hadronic uncertainties. 
Studying their low q 2 dependence can be helpful in discerning among various possible 
new physics scenarios. All three asymmetries can be extracted from the full angular 
analysis of B — > K*£ + £~ . Our formulas apply to both the massless and the massive 
lepton case. 



PACS: 12.20.He 



1 Laboratoire de Physique Theorique est une unite mixte de recherche du CNRS, UMR 8627. 



1 Introduction and basic formulas 



Searches for physics beyond the Standard Model (BSM) through the low energy flavor 
physics experiments could provide us with an indirect tool for identifying the new particles 
that will hopefully be directly detected at LHC. Flavor changing neutral currents are 
known to be particularly revealing in that respect, and the b — > s transitions are known 
to be particularly interesting. In this paper we will focus on the exclusive b — > s mode, 
B — > K* (Kir)£ + £~ , where £ stands for a lepton. Q After an intensive research devoted to this 
decay it has been understood that the full branching fraction is extremely difficult to handle 
theoretically because: (1) all the hadronic form factors enter the corresponding expression, 
and (2) every form factor should be integrated over a very large range of kinematically 
accessible values of q 2 , namely q 2 G [4m 2 , g^J, with g^ ax = {m B — m K *) 2 = 19.25 GeV. 
None of the available methods to compute form factors is viable in the entire physical range 
of momenta transferred to leptons, and therefore only the partial decay rates (integrated 
over a relatively short span of g 2, s) can be computed, albeit with a limited control over 
theoretical non-perturbative QCD uncertainties. Furthermore, since the study of this decay 
is based on the use of operator product expansion (OPE), it is essential to avoid the 
resonances arising when the energy of the lepton pair £ + £~ hits the production thresholds 
of the cc- resonances (i.e. J/ip and its radial excitations). For that reason one should try 
and work below q 2 = m 2 ,^ ~ 9.6 GeV 2 , or experimentally veto the points at which the 
resonances are expected to occur. Besides, one can also encounter difficulties with the 
charmless resonances (p, u> and their radial excitations) but that effect turns out to be 
practically negligible thanks to the CKM suppression. 

Using OPE, at low energies, this decay is described by the following effective Hamilto- 
nian pQ: 



%cff — 7^~VtbVts 



.i=l i=7,8,9,10,P,5 ^ ^ 



where the twice Cabibbo suppressed contribution (oc V u bV* s ) has been neglected. The 
operator basis looks as follows [21 13] : 

1 = (s 7M T a P L c)(c7 M T a P L 6), 

2 = (s^P L c)(c^P L b), 

q 

O, = (rs^T a b)J2(^T a q), 

i 

practice one considers either electron or muon, but not tau, although the expressions given in this 
paper apply equally to the B — > K*t + t~ case. 
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O r = -^{m^P^F^, 0' 7 = - 2 m b {m, v P L b)F^ , 

g 2 g 2 

8 = -m b {sa, u T a P R b)G^\ 0' 8 = -m b (sa^T a P L b)G^ a , 

g g 

2 2 

9 = ^(s^P L b)(h»£), <D' 9 = e -{s llx P R b){hn), 

2 2 

Oio = - 2 (s^P L b) (hV) , O' 10 = - 2 (s^P R b) (Ind) , 

g 2 g 2 

2 2 

O s = -^—m b (sP R b)(££), 0' s = -^ mb (sP L b){i£), 

1071" l07T^ 
2 2 

Op = ^m b (sP R b)(£ l5 £), 0' P = ^m b (sP L b){£ l5 £), (2) 

with P^ p = (1 =f 75)/2. The operators with the opposite chirality, 0[_ 6 , have been 
neglected. Short distance physics effects, encoded in the Wilson coefficients Cj(/i), have 
been computed in the Standard Model (SM) through a perturbative matching between 
the effective and full theories at n — mw, and then evolved down to the \i — m b by 
means of the QCD renormalization group equations at next-to-next-to-leading logarithmic 
approximation (NNLO) [2]. For the reader's convenience the resulting values of Wilson 
coefficients Ci(m b ) are listed in Appendix of the present paper. An important feature 
regarding the operator Og is its mixing with Oi,...^ through diagrams with a virtual photon 
decaying to £ + £~ . It is customary to reassemble Wilson coefficients multiplying the same 
hadronic matrix element into effective coefficients appearing in the physical amplitudes, 
namely jl] 

rcS 47T 1 4 20 80 
47T 1 10 

Cg = — Cg + C3 — - C4 + 2OC5 — — C/6 , 
a s 6 3 

Of = -C 9 + YV), 

^yeff _ 4-7T / cff _ 4tT ^, . . 

°10 — a °10 5 °7,8,9,10 ~~ a °7,8,9,10 ' 



where the function Y(q ) is 



4 64 64 1 / 4 64 

= 3C3 + y^ + ^C 6 -2% 2 ; 0)^3 + -C 4 + 16C 5 + y a 

/4 \ 1 / 4 64 

+% 2 ,m c ) f-d + Ca + eCa + eOCsJ --% 2 ,m fe ) ( 7C 3 + - C 4 + 76C 5 + — C £ 



(4) 



and 



, 1 arctan ; 2 > 1 

4 / m 2 2 \ 4 . ^/z 



(ui 2 \ 4 
In - 3 - z\ - - (2 + z)y/\^\\ x 



9\Ai 2 3 / 9 v ;VI 1,1 + VI - z i7T 

m ^ z < 1 

2 

(5) 
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with z = Am 2 /q . As far as the long distance physics effects are concerned, they are 
encoded in two hadronic matrix elements that are conveniently expressed in terms of seven 
Lorentz invariant form factors. The hadronic matrix element of the standard V-A current 
is decomposed as, 



(K*(k)\s ltl (l- l5 )b\B(p)) = e^ pa e* v pPk° 2V{ - q > - ie*Jm B + m K .)A x {^ 

m B + m K * 



M{p + k)»{e* ■ q) ^ + iq»(e* ■ q) ^ [A 3 (q 2 ) - A (q 2 )} , (6) 
ttl b -\- nix* q 

where is the .^-polarization vector, and thanks to the partial conservation of the axial 
current, 

a / 2\ m B + m K * 2 m B - m K * 2 

also ensuring that no artificial divergence emerges at q 2 = 0. Other 3 relevant form factors 
parameterize the matrix element of the electromagnetic penguin operator, 

{K*{k)\8<r^{l + l6 )b\B{p)) = 2* W e*V^T 1 (g 2 ) 

+ [eliml - m\») - (e* ■ q) (2p - q) J T 2 (q 2 ) 



+ ■ q) 



q 2 

% 2 T~ (P + k )^ 

m B — m K * 



Ts(q 2 ), (8) 



with Ti(0) = ^2(0), ensuring that only on form factor describes the physical B — y K*j 
decay. 

The explicit expression for the full differential B — y K*{— y Kit)£ + £~ decay rate was 
presented in ref. [U] triggered a quite intense activity in searching for good observables, 
namely those that simultaneously satisfy three requirements: (i) to not suffer from large 
hadronic uncertainties, (ii) to have a pronounced sensitivity to the presence of physics 
BSM, and (iii) to be experimentally accessible at LHCb and/or Super-B factories. Written 
in terms of four kinematic variables, the differential decay rate reads [6] (see also [3]): 

d 4 T{B° K*°£+f) 9 2 



dq 2 d cos 9i d cos 9 k d<p 32n 



I(q*,9 e ,0 K ,<f>), (9) 



where 



/(g 2 , 9 t , 9 K , 0) = I{ (q 2 ) sin 2 9 K + I{{q 2 ) cos 2 9 K + [l s 2 {q 2 ) sin 2 9 K + / 2 V) cos 2 9 K ] cos 29 t 
+ h(q 2 ) sin 2 6 K sin 2 9 e cos 20 + I^q 2 ) sin 29 K sin 29 £ cos 
+ J 5 (g 2 )sin2^ sin ^ cos (10) 
+ [/g(g 2 ) sin 2 6^ + /^(g 2 ) cos 2 6^] cos 6^ + I 7 (q 2 ) sin 20^ sin ^ sin 
+ Is(q 2 ) sin 26*^ sin 26^ sin0 + Ig(q 2 ) sin 2 0# sin 2 0£ sin 20 . 

Besides q 2 , the other kinematical variables are defined with respect to the direction of flight 
of the outgoing K°* in the rest frame of B°. In particular 9% is the angle between that axis 
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and lepton i~ in the rest frame. 9k, instead, is the angle between that same axis and 
K~ in the Ktt rest frame. Finally (f) is the angle between the Ktx and £ + £~ planes (see 
ref. ??). PI The functions Ii(q 2 ) are related to the transversity amplitudes >lj_,|[,o,s,t(<Z 2 ) as 
follows: 



where 



4m« 



KI + Wl + q . 



[\A t \ 2 + 2Re(A%AK*)]+ff\As\ 



ft 
4 



.4 



L |2 



1-4 



L|2 



1 



-$(K L I 2 + IAT) , 
/3l[|Ai| 2 -|4i 2 + (L^it:)] , 



/3| [Re (4^f ) + (L -> R 



V2 



V2p t 



Re (4^1* -(L^-R)) 



Re (Af A* 5 + AfA* 8 ) 



2f3 e [Re (Af4*) -(L->R 



[Re (4^) + (L -+ R)] 



1 

71 



Im (44* - [L -> i?)) + 
/3| [Im (44*) + (L ->• R 



Im (AiAJ + A*A* S ) 



ft [Im (44*) + (L -> i? 



4m? 



(11) 



(12) 



while the amplitudes A 



,0,SU 



.4 



0St (g 2 ), when written in terms of the Wilson coeffi- 



2 One should be careful and distinguish the CP conjugated mode B° — > > K + ir~)£ + £~ . 9i and <fi 

are defined in the same way and are related to the ones discussed in the B° — > K°*(— > K~ir + )£ + £~ via 
9{ — > 7r — 9e, <f> — > —<j). The net effect on the coefficient functions in the angular distribution <fToj) relevant 
to our discussion in this paper is that Ie(q 2 ) and Ig(q 2 ) would change the sign. See refs. [51 117j. 
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cients and the form factors, read: 



2m 



b (Cf + Cf) 7\(g 2 



1 



+ [(Cf + Cf) =F (C£ + C 
2m b 



eff/ 
10 



V(q 2 



m B - m K , 



vtib + ttik* 

cS'\ rr I 2\ 



(cf -cf')7 1 2(g 



+[(cf -cfOTtc'i'o-C': 



eff/ 
10 



^(g 2 



tfi5 - m^* 



(13) 



(14) 



4<*(^ 



A(g 2 ) 

A s {q 2 ) 



2m,K* v <? 2 



x 



(m| - m^, - q 2 )(m B + m K *) Ai(q 2 ) - X(q 2 ) 



A 2 {q 2 



+2m b (Cf - Cf) 



(m 2 B + 3m^, - q 2 ) T 2 (q 2 



vtiB + rriK* 
A(9 2 ) 



m B — m 



-T 3 (q 2 



K* 



IAf\(q 2 



2 (Cf - C[f) + ^-(C P - C' P 

mg 



A (q 2 



-2y^f\(q^(C s -C' s )A (q 
In the above formulas, q 

Af=\V tb V*\ 2 



3 2 10 7r 5 m| 



q 



KQ ) = [Q - ( m B + m K *) } [q - (m B - m K 



(15) 
(16) 



(17) 
(18) 



(19) 



2 Transverse asymmetries 

After integrating eq. ([£]) over Qg and 9k, one arrives at: 



d 2 T{B -> K*£ + £~) _ 1 dT 



1 + l -F T (q 2 ) (A {2 \q 2 ) cos 20 + A^\q 2 ) sin 20 



<ig 2 (i0 

where Fplq 2 ) is a fraction of the decay product with transversely polarized K* 

F(n2) # l^(g 2 )| 2 + |A||(g 2 )| 2 



(20) 



(21) 



3 We chose to work in the basis of operators given in ref. [3] in which Os,p and 0' s P are dimension 7 
operators, and therefore the corresponding Wilson coefficients are not dimensionless but of dimension — 1, 
and therefore the two terms within the brackets in A t (q 2 ) are of the same dimension. 
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2.1 A%\q 2 ) 

The transverse asymmetry A^\q 2 ) in eq. (20), first introduced in refs. [SI Ej, is CP- 
conserving and it reads 

where \A\\ ,_i_(g 2 )| 2 = ^^(g 2 )] 2 + ^^(g 2 )) 2 . That quantity is expected to be experi- 
mentally accessible at LFICb and Super-B factories, and nowadays is considered as one 
of the most interesting observables to study since it satisfies all three above-mentioned 
requirements. What is interesting to note is that A^\q 2 ) only involves A^j_(q 2 ), and not 
A),s,t(g 2 ) amplitudes. The latter ones are much more difficult to handle in QCD. In many 
phenomenological studies the symmetry relations among form factors, first demonstrated 
in ref. [8], are used to express all the form factors in terms of only two functions, £j_(g 2 ) 
and £|| (g 2 ). The problem with that approach is that this approximation is only valid in the 
limit of m B — > oo (i.e. E K * — > oo, at low g 2 's). It is not clear how to compute the non- 
perturbative 1 / corrections to that approximation, as well as those due to the finiteness 
of rriK*- Perturbative corrections can be handled in the QCD factorization approach [9] 
or via the soft collinear effective theory [10], but we do not know how to compute the 
non-perturbative corrections from principles of QCD. In particular it is not clear whether 
or not it is even possible to formulate the problem and compute the relevant 3-point cor- 
relation functions through numerical simulations of QCD on the lattice, and then extract 
the needed £y(g 2 ) and £_i_(g 2 ) form factors. Instead, the full form factors, as defined in 
eqs. (6j8 ), can and have been computed on the lattice. Particularly difficult appear to be 



the form factors T 3 (g 2 ), A 2 (g 2 ), and A (q 2 ), that are linear combinations of £j_(g 2 ) and 
£|| (g 2 ) as soon as the mass of K* is not neglected [see eqs. (107,109,113) in ref. [8]]. On 
the lattice the signals for these three form factors are particularly difficult to keep under 
control [see for example v4 2i0 (g 2 ) in refs. [TH [12]]. The advantage of using the quantities 
involving only Ah t ±(q 2 ) is that they do not require a detailed knowledge of these three form 
factors. Moreover, the ratios A 1 (q 2 )/T 2 (q 2 ) and V A (g 2 )/T 1 (g 2 ) seem to verify the symmetry 
relations of ref. [SJ and exhibit a flat g 2 -dependence over a wide range of g 2 's, satisfying 
the approximate relation |^] 

A 1 /T 2 _ V/Ti ^ Q2GeV -i (23) 



rag — nix* rriB + tuk* 



4 Since these ratios are flat for low <? 2 's, from now on we will keep the argument of these form factors 



implicit. In estimating the ratios (23) we take V/T\ ~ 1.24 and A\jT-z ~ 0.88 from ref. [TJ], which is 
consistent with what is obtained through an alternative QCD sum rule method, i.e. V/T% ~ 1.22 and 
A1/T2 — 0.92 [TJ], although the results of these two methods for the absolute values of the form factors 
are quite different. One of us (D.B.) checked that this feature is also verified in (quenched) lattice QCD. 
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Finally when replacing ^4ii_|_(g 2 ) in eq. (22) it is convenient to factor out T\(q 2 ) which then 



cancels out in the ratio. Schematically that amounts to 



4 2 V) 



V/T x 



m B + m K * 



en + a 



A x jT 2 



m B - m K * 



C ± + C\ 



V/T, 



m B + m K * 



Cw + a 



AJT 2 



m B - m K » 



(24) 



where, for notational brevity, the combinations of Wilson coefficients and kinematical fac- 
tors are denoted as C^,,. This form is useful because the ratio T 2 (g 2 )/Ti(g 2 ) is a well 
controlled quantity and is 1 at q 2 = by definition, i.e. 



T 2 (q 2 ) T 2 (0) 2 2 

xW) = W) +zq =1 + zq 



(25) 



where 2 is a slope which in a simple pole model is given by z ~ — l/m^, = —0.034 GeV 2 . 
That value is close to z = —0.028 GeV -2 , as inferred from the light cone QCD sum 
rules [13J , as well as to the one obtained in (quenched) lattice QCD, z = -0.030(3) GeV~ 2 [15]. 

Therefore the advantage of using the quantities involving the amplitudes A^±(q 2 ), and 
not AQs,t{(f)) is that the relevant hadronic uncertainties are under better control. The 
price to pay when dealing only with ^4||,j_ (^ 2 ) is that no information about a coupling 
to the new physics scalar sector can be accessed. This, of course, can be viewed as an 
advantage too, because the number of possible new physics parameters to study becomes 
smaller. A particular interest in studying A^\q 2 ) is the fact that its value at q 2 = is 
identically zero in the SM, while in the extensions of SM in which the couplings to right 
handed currents are allowed (oc C' 7 ) it can change drastically, 



4 2) (o) 



2 Re [CfC'f*] 
\Cf\ 2 + \C' 7 eS \ 2 



(26) 



This argument is strictly valid only at q 2 = 0, while away from that point the situation 
becomes more complicated since A^\q 2 ) receives non-negligible contributions from the 
terms proportional to Cg^ and Cg e f . 

Written in terms of 7j(g 2 )'s explicitly given in eqs. (9 10 11), A^\q 2 ) becomes particu- 
larly simple, 



4V 



2/f (g 2 



(27) 



It is important to note that there is an identity I{{q 2 ) = Sl^q 2 ), valid in the case of massless 
leptons (reasonable assumption for B — > K*e + e~). Lepton mass breaks this identity but it 
remains a good approximation for q 2 Am 2 (which in practice means a couple of GeV 2 ). 
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2.2 



4 m V 



Measuring the term proportional to sin: 
quantity 



in eq. (20) is highly interesting because the 



2 Im 



AUq 2 )A L ± *{q 2 ) + A*{q 2 )Af ( 



\A ± (q 2 )\ 2 + \A {l (q 2 )\ 2 



(28) 



is obviously identically zero in the SM for all accessible q 2, s. It may acquire a non-zero 
value only if there is a new phase, coming from physics BSM, and therefore its non-zero 
measurement would be a new physics discovery. In particular, at q 2 — > one obtains 



4 m) (o) 



2 Im [CfC'f*] 
\Cf | 2 + \Cf s \ 2 



(29) 



so that its non-zero measurement would mean that the new phase comes from the elec- 
tromagnetic penguin operator, either C-j or Cj, or both (if they are different). As in the 
case of A^\q 2 ), the above formula at q 2 = should be taken as asymptotic, because away 
from that kinematics other Wilson coefficients {Cg S 10 and Cg c f ) become important too. 

From the g 2 -shape of A ( '^ a \q 2 ) one can get some valuable information about the source of 
the new physics phase, as we will discuss in the next section of the present paper. No- 
tice also that our definition of A^ m \q 2 ) is different from Ai m (q 2 ) introduced in ref. |16j : 
2fJ 2 AUq 2 )=A^\q 2 )F T (q 2 ). 

A^°\q 2 ) too has a particularly simple form when written in terms of angular coefficient 



functions defined in eqs. d9fl~0fTTj ), namely 



4 m V) 



2/|(g 2 



(30) 



2.3 



4\ q 



2 " and Ap\q 2 ] 



Another quantity that involves only the amplitudes Ah t ±(q 2 ) has been recently proposed 
in ref. [TBI: 



4V 



\A[Aq 2 )Af(q 2 ) + Ai(q 2 )A**(q 2 ) 



\A±(q 2 )\ 2 + \AAq 2 ) 



(31) 



It was then shown in ref. [T7] that A^\q 2 ) arises naturally after realizing that eq. (9), 
in the limit of massless leptons (me = 0), respects 4 symmetries. More specifically, the 
coefficient functions h(q 2 ) in (10) are invariant under the following transformations: 



o Two independent phase transformations, 



A V 



Af(q 2 ), A?\q 2 ) = e^A?{q 2 ) (*=±,||,0) 



(32) 



8 



o Two rotations of suitably defined vectors, n\ = (A^,A^*), n 2 = {Aj_,— A^*), n.3 
(A%,A**), 



R(9) 



cos 9 — sin 9 
sin 9 cos 9 



V{9) 



cos 9 —i sin 9 
-i sin 9 cos 9 



(33) 



They eventually showed that for the massless lepton, ^ 

l( %o, Vie/I V) - WW) ~ 4(V(g 2 ) + h\q 2 )) 



8/|(g 2 



(34) 



m,£=0 



Since we consider only the quantities constructed from A^'^(q 2 ), or better the observables 
involving I 26 (q 2 ) and I 3t g(q 2 ), it suffices to consider only the first two vectors, n\ and n 2 , 
in which case the reasoning used in ref. [T7j can be extended to the massive lepton case. 
Using the two global phases, 4>l and (f>u, one can first assure that A^' R (q 2 ) are real and 
positive. Furthermore, by choosing a suitable angle 9, one can rotate A^ away and have 
ni = (0, A R ). The scalar products of vectors ni^ can then be easily expressed in terms of 
functions 7j(g 2 )'s as follows: 

o 2/|(g 2 ) - I 3 (q 2 ) 



K| 2 = (Af (g 2 ))- 



n 2 



2 = \A L ± (q 2 )\ 2 + \A R (q 2 )\ 2 = - l(g2) + hiq2 



ni -n 2 = -A R {q 2 )Af{q 2 ) 



PI 

f3dl(q 2 ) + 2ih{q 2 



(35) 



or equivalently, 



A\{q 2 ) 



Af(q 2 ) = 
Af(q 2 ) = 

\A L ± (q 2 )\ 2 



2P 2 {q 2 ) - I 3 (q 2 ) 



1 (3 t Il(q 2 ) - 2il 9 (q 2 



1_ 16/|(g 2 ) 2 -4/ 3 (g 2 ) 2 -/3 2 / 6 s (g 2 ) 2 



4/ 9 (g 



2\2 



4/3| 



2II(q 2 ) - I 3 (q 2 ) 



(36) 



After inserting the last expression into eq. (31 ) we then simply obtain 



A^(q 



iW|V) - p 2 M 2 )i?(q 2 ) - 4(/ 3 V) + / 9 V)) 



(37) 



8/|(g 2 ) 

This formula is valid for any mass of the lepton, and obviously reproduces the massless 



lepton formula (34) derived in ref. [T7]. One should be very careful in constructing the 
observables from Ii(q 2 ), and not include those that do not respect the above symmetries. 



'Notice that in the massless case If{q 2 ) = 3/f (g 2 ), which is not valid if mi ^ 0. 



9 



The observable A^\q 



2 ^ is indeed independent from A^\q 2 ) 



and 4 m) (g 2 



In fact it 



is easy to see that there can be at most three independent asymmetries built up from 
A±i\(q 2 ), respecting the above symmetries. In a generic scenario of physics BSM there can 
be 8 real amplitudes: A^ R (q 2 ) and A^' R (q 2 ), each with its real and imaginary parts. The 
above four symmetries help reducing the number of independent combinations to four. |^] 
Since we consider the ratios, that implies the restriction to at most three such independent 



observables. In our case those three are: A^\q 2 ) 



From eq. (37) it is clear that A^\q 2 ) is too complicated a quantity. After introducing 



4 m V) 



and AP(q 2 ) 



A { 4\q 2 ) and A ( ™\q 2 ), the only new angular coefficient function needed for A T (q 2 ) is 
Jg(g 2 ). We find therefore more convenient to introduce, 



4 e V 



2 Re 



Ak{q 2 )Al*(q 2 ) - A R {q 2 )Af{q 2 



|A ± (^)|2 + | A|| (g2)|2 4/|(g2) 
4 V) is related to A {2 \q 2 ), A^ m) (q 2 ) and A { ™\q 2 ) via 

[2AP(q 2 )) 2 + (4 2 V)) 2 + (4 m V)) 2 + (4 re V) " 



(3f 



1 . 



(39) 



It is worth mentioning that in the scenarios in which the new physics does not modify the 
scalar amplitude, A s (q 2 ) [c.f. eq. (18)], the angular coefficient function I§(q 2 ) = 0, and 



At(q 2 ) is related to the usual forward-backward asymmetry as 



4 re V 



F T (q 2 



(40) 



where Fxiq 2 ) = 4/|(g 2 ) / (dT / dq 2 ) is the quantity we already defined in eq. (21), whereas 

d 2 T r0 

I d cos t>£-j—. 

A FB (q 2 ) 



f dcosOt-— - 

o dq z d cos 6v 



d 2 Y 

dcos9l T^ Z~ 

_] dq z d cos Vi 



3 2/ 6 V) + / 6 c (^ 



dT/dq l 



dT/dq 2 



(41) 



The expected shape of A FF ,(q 2 ) in the SM and its various extensions has been discussed 
in great detail in the literature [T5]. Importantly, the property that limA FB (g 2 ) — > 0, 

<j 2 -s>0 

caries over to A^ e \q 2 ). Therefore, unlike A^ (0) and 4™ (0) whose values at q 2 = 
can change considerably if the new physics affects the Wilson coefficients C\ , the third 
quantity Aj, e \o) remains insensitive to new physics. However the g 2 -shapes of our three 
asymmetries can teach us about new physics. In particular the point at which these three 
quantities might become zero can be helpful in discerning among BSM scenarios. 

2.4 How to extract A^ e \q 2 ) from experiment? 

Before passing onto the potentially interesting phenomenology that one can learn from 
the above asymmetries at low g 2 's, we need to discuss the possibility to experimentally 



3 In terms of angular coefficient functions, those four functions are: I|(<? 2 ), I^iq 2 ), I|(g 2 ), and Ig(q 2 ). 
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measure the above introduced asymmetry A^ c ' 



.,</ ) [and therefore A^\q 2 ) too]. Since 



A^ e \q 2 ) resembles the standard forward-backward asymmetry, a way to extract it from 
experiment is quite similar to that employed for Afb(<? 2 )- After integrating eq. ^ in 
one can separate the events with the lepton going forward from those in which lepton flies 
backwards and get the differential distribution of the forward-backward asymmetry, 



d 2 A FB {q 2 ,9 K ) 
dq 2 d cos 6 k 



d cos 9t 



2tt 



d 4 T{B -> K*i 



dq 2 d cos 0£ d cos 9 k d(f) 
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\Q{dT/dq 2 



F T {q 2 



9 



16/3, 



dV/dq 2 

[l°(q 2 ) sin 2 9 K + I c 6 (q 2 ) cos 2 9 K ] 



A^(q 2 )sm 2 9 K + fr 



4J|(g2 



(42) 



As before, -Fr(<? 2 ) is the quantity specified in eq. (Elk, and Aj e) (g 2 ) can be simply read off 
from the angular dependence of the forward-backward asymmetry. The term proportional 
to cos 2 9k is expected to be small: in fact in the SM I§(q 2 ) is exactly zero, while in the BSM 
scenarios it can become important only if the new scalar contributions are large. Q However, 
even in that latter case I§{q 2 ) remains negligibly small if one considers B — > K*e + e~, 
because the scalar amplitude As(q 2 ) enters in I^q 2 ) multiplied by the lepton mass. In 
the B — > K*ji + fi~ mode, instead, one should check whether or not the term cx cos 2 9k is 
discernible. 



3 What can we learn from the shapes of Aj, (q 2 ) 
and from A^ 2) ' (re) ' (im) ( % 2 ) -> 0? 

As we already mentioned, the quantities A^ m \q 2 ) and A^\q 2 ) can considerably change 
the value at q 2 = in the presence of new physics contributions. A^ e \q 2 ) instead scales 
as q 2 at low g 2 's, and is zero at q 2 = in the SM and its extensions. Even the functional 
dependence of A^ e \q 2 ) is not too sensitive to the presence of new physics, which is different 
from what happens with A^ a \q 2 ) and A^\q 2 ). It is interesting to examine the situations 
in which these asymmetries can change the sign, while remaining in the low g 2 -region. 

To measure a non-zero value for the asymmetry A^, (g 2 ), a new phase is needed. The 
origin of that phase can be figured out from the functional dependence of A t ^ ra \q 2 ) at low 
g 2 's. In particular A^ a \q 2 ) can go through zero at q 2 . ^ if there is a solution to the 

7 From the experimental upper bound B(B S — > [121 122] one can constrain \Cs\ < 0.8/to& < 

0.2 GeV" 1 . 
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Figure 1: Transverse asymmetries in the Standard Model (C' 7 = C' 9 = C' w = 0). Note that A^ m \q 2 ) 
(dotted line) is identically zero, A^\q 2 ) (full line) indistinguishable from zero for low values of q 2 's, while 
A^(q 2 ) varies between 1 and — 1. 



equation: 



Im <{ Cf'Cf* - 



1 q 2 



+ 



2mi 



2 2m h 



(Cf - Cf')(Cf + ci 
+{Cf + Cf)*{Cf 



m B + m K * 
9 m B - m K , 



°9 "T °10 °10 



V(q 2 )lT l (q 2 )A 1 (q 2 )/T 2 (q 2 
triB + rriK* tub — tuk* 



. (43) 



On the other hand, measuring the point at which A^, e \q 2 ) crosses the g 2 -axis resembles 
the standard discussion of A FB (ql) = 0, and the corresponding q^ is defined via 



nQ 1 V°10°9 — °10 U 9 ) 



Ai(g 2 )/T 2 (g 2 ) VV)/Ti(g 2 ) 
ttib — itlk* rriB + rrix* 



ICS 



+ (^7 

+ (cf +cf )*(c; 



C cff/ )(C cfr + (7 cfr/ )* ^"(^V^iX? 



eff 
10 







A 1 (q 2 )/T 2 (q 2 ) 
m B - m K * 



0. 



A similar formula for determining q^ , 2 s at which A 1 ^' (q 2 



(44) 



' ' 1 ' crosses zero cannot be written 
in a simple form. To find a simple formula, it suffices to note that for small g 2 's a following 



approximation is reasonable, Ti(q 2 )\ l l 2 (q 2 ) m T 2 (q 2 )(m 



B 



m 



K* 



), so that A^p 



(q 2 



can 



12 



be written as 



r Re 



mr 



{Cf + Cf){Cf + Cf) 



rriB + rriK* 



(Cf-Cf')*(C c 9 



m B - m K , 



2m 2 



A 1 {g 2 )/T 2 {g 2 

m B - m K * 



2 



[\cf+cf\ 2 + \c$- 



°io I J 



\Cf 



+ \c 



10 



^cfr/121 
°io I J 



+ 4Re(CfCf* N 



0. 



(45) 



Furthermore it is reasonable to consider Cg ff (g 2 ) to be g 2 -independent at low g 2 's, because 
in that region the function Y(q 2 ) in eq. (|4|) varies very slowly at low g 2 -region. In fact, 



Y(q 2 ) makes about 20% of the whole Cf\q 2 ), so that its small variation in q 2 might at 
worst entail an error of a couple of percents. In the following we will therefore consider 
Y(q 2 ) = 1. To situate the positions of g 2 q 2 ^ and q 2 , 2 n it is helpful to use eq. (23), 
that we will denote by R. The expressions for q 2 ( 43[44|45 ) simplify, but their solutions are 
not very compact unless we restrain our attention to the specific NP scenarios. We should 
emphasize once again that the discussion above refers to the case in which C 7 remains 
negative, i.e. of the same sign as in the SM. If that was not the case, A J Te \q 2 ) would have 
not crossed the g 2 -axis. 



3.1 New physics in CV and C 7 

From now on, for notational simplicity, we will drop the superscript "eff" from the Wilson 
coefficients. If C$ and Cio are kept fixed to their SM values, and C' g = C' w = 0, the 
points at which the asymmetries A J 2 \q 2 ), A t >^ a \q 2 ) and A J T , e \q 2 ) may change the sign are 
obtained by solving eqs. ( 45p3 , 44) respectively. We obtain: 



<?0 (2) 



2m b Re{C;C 7 ) 
' R C 9 Re{C 7 ) 



2 2m b lm(C^C 7 ) 

<?o(im) - - R CghaJOr) ' 

2 2m fe Re(C 7 ) 

%M = -^R~cr— ■ (46) 

where CV, C 7 G C is understood. At this point we should reiterate that the approximations 
discussed above |^] are used to obtain the simple expressions for q 2 . q 2 ^ and q 2 but 
all the plots presented in this work are obtained by using full expressions. Excellent agree- 
ment between the approximate and full results confirms the validity of our approximations. 
One can distinguish three interesting situations: 



J For reader's commodity we repeat that our approximations at low q 2, s are: (l)(Ai/T2)/(rriB — mx* 
T 1 )/(m B + m K ,) = R, (2) dCf/ 
T 1 (g 2 )A 1 /2( g 2 ) = T 2 (q 2 )(m% - m? K .). 



(yfTx)f{m B + m K *) = R, (2) dCf/dq 2 = 0, and Y(q 2 ) = 1, and (3) in the case of {q 2 ) we also use 
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(c^\-\Cf/2\e i7T ^ 





^-|Cf||e^/ 6 ,Cf I /2^ 




Figure 2: Dependence of the asymmetries A^\q 2 ) (full line), A { ™ l] ' (q 2 ) (dotted line), A^ (q 2 ) (dashed 
line) on q 2 for different possibilities (real or complex) of (C 7 ,C 7 ), while (Cg,Cg) and (C\q,C' w ) are kept 
fixed to their SM values. 



For C 7 , C 7 G E, and C' 7 ^ 0, one has 

1o (2) = Qo (rc) = ~"]=r^r , 1o (im) - does not exist . (47) 



If C 7 = 0, the situation becomes the SM-like, namely ,s remains the same, while 
/ 2 x does not exist anymore (see also fig. Fij). 



Particularly interesting is the situation with C 7 G M and C 7 G C, in which all three 
asymmetries cross the g 2 -axis at the same point, 

2 2 2 2m;, C 7 . . 

% (2) — % (im) — % (re) ~ JT Q) ' 

Of course this coincidence of zeroes could be spoiled if our approximations were bad. 
We checked that three asymmetries indeed become zero at almost the same point 
even when the full expressions are used, as shown in fig. [2] 
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• In a reverse situation, i.e. C 7 G C and C' 7 G R, one has 

2 2 2mf,Re(CV) 21 • / A ^\ 

% (2) = % (re) = ^ (T— , % (im) ~ d ° eS not eXlst ■ ( 49 ) 

It is important to note that the above discussion holds in the situations in which the new 
physics does not alter the sign of the Wilson coefficient C7. If that happens then obviously 
the asymmetries do not have zeroes anymore. F] We stress again that in order to produce 
the plots in fig. [2] we did not use any approximation to the full formula. To illustrate the 
situations in which the new phase alter the SM values of the Wilson coefficients, we choose 
that phase to be 7r/6, but without altering the sign of Re(C 7 /) ) with respect to the SM 
value. In other words, we take for C 7 G C, C 7 = -\C* M \e in/6 = Cf M e i7r/6 . Similarly, when 
C' 7 G C, the illustrations are provided by using C 7 = -|Cf M |/2e W6 . The bands of the 
curves shown in all the plots of this paper are obtained from Monte Carlo in which we 
uniformly varied the form factor ratios R, the slope z, and the quark masses within the 
ranges specified in Appendix. 

3.2 New physics in Cg and Cg 

Another distinct possibility is to consider the case in which Cg and C'g are affected by 
the new physics effects, while C 7 and C w remain intact (and C' 7 = C' w = 0). The points 
#0 (2) (im) (re)' wnere the asymmetries have zeroes, are: 



^0 (2) 



2m b C 7 Re (C' 9/ 
"~R Re (C 9 C 9 ) 



2 2m b C 7 Im (C 9 ) 

% (im) — 



R Im (CgC'g, 



2 2m fe C 7 

% ^ ~ R~RMCgj' (50) 

Three interesting situations in this case are: 

• If there is no new phase, i.e. Cg, C'g G K, one has 

?o (2) = ll (rc) = ""^^ ' ^0 (im) - does not exist . (51) 



For Cg G R and C'g G C, the three zeroes again coincide, 

2 2 2 C-j 

% (2) — % (im) — % (re) ~ TT 7T • l 0Z J 



9 That issue has not been resolved by the recent results on Ayb{<1 2 ) measured at BaBar and Belle fZT\ . 
10 At q 2 > 4m 2 ., Im[y(g 2 )] ^ 0, and therefore Cg becomes complex. That region is however above the 

n,2 



point at which the asymmetries cross the q -axis. 
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-1.0 



Figure 3: Similar to 6g-\^ plotted is the q 2 -dependence of three asymmetries A^ (q 2 ) (full line), (q 2 ) 
(dotted line), A^(q 2 ) (dashed line), for different possibilities of new physics in (Cg,Cg) (see the text), 
while (CV, C' 7 ) and (Cio, C 10 ) are kept fixed at their SM values. 



• If, instead, Cg G C and C 9 G K, then 

2 2 2rrib Cj 2 . . 

% (2) — % (re) — ^~R e (C 9 ) ' ^ (im) ~~ ' ^ ' 

These three situations are illustrated in fig. [3j We see that, like in fig. |2j when one of 
the two Wilson coefficients is complex while the other one is real, two or three asymmetries 
cross the g 2 -axis at the same point. That information alone would not tell whether the new 
physics contribution modifies (ffl or Cg . However the low g 2 's shapes of the asymmetries 
A%\q 2 ) and A^ m \q 2 ) in fig. E are very different from those in fig. 2 which solves the 
ambiguity. 

3.3 New physics in C\$ and C' 1Q 

The third special case we consider here is the one in which C\$ and C' w receive non-zero 
contributions from new physics, while CV and Cg retain their SM values. The zeroes of our 
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Figure 4: Expected q 2 -shapes of asymmetries A^\q 2 ) (Ml line), A^ m '(q 2 ) (dotted line), A^ e \q 2 ) (dashed 
line) , if the new physics alters only (Cio, C' w ), and (C7, C' 7 ) and (Cg, C 9 ) remain at their SM values. 



three asymmetries then become: 

2 2 _ n 2 

% (2) — % (im) — U > % (re) 



2m;, C7 
i? Cg 



(54) 



independent on Cio and C' w . 

In this case, neither A^\q 2 ) nor v4^ m ^(g 2 ) crosses the g 2 -axis away from zero. The q 2 
dependencies of the three asymmetries are different from those found in the previous two 
cases, as it can be appreciated by comparing fig. [4] with those presented in the previous 
two subsections. 

3.4 Maximum of Ap\q 2 ) 

In all this discussion we supposed the new physics would not change the sign of the real 
part of the Wilson coefficients 6*7,9, i-e. their signs in the SM. From the plots presented 
in figs. TJ 2 3 and 4, we see that the shape of A^ c \q 2 ) remains stable under the variation 
of Wilson coefficients, even in the presence of a new physics phase. Its maximal value 
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can provide us with another interesting information. Using the approximations discussed 
above, in the SM, we get that the point q 2 at which dA^ e \q 2 ) / dq 2 = is 

= ^ - 1-3 GeV 2 4V) = 1. (55) 

We checked numerically that, even without using our approximations, this result remains 
valid. In the presence of new physics, the value of q 2 only slightly shifts with respect to the 

SM one, but 4 re V) 

becomes lower. The only exception to this rule is the situation in 
which the right-handed currents are absent and all the Wilson coefficients receive the same 
phase. Concerning the three new physics scenarios discussed in this section, we note that 
Aj, e \cp) is considerably lower in the case of new physics modifying the (C7, C' 7 ) values: for 
(C 7 , Cj) = (C 7 SM , e^|C 7 SM /2|), one gets A^ffi) w 0.6 regardless of the value of the phase 
ip; for (C7, C 7 ) = (— |Cf M |e^, C 7 M /2), the suppression of this asymmetry is even stronger 
but depends on (p [in particular, for tp = rr/6 we have A^icj 2 ) ~ 0.4]. 



4 Summary 

In this paper we discussed the benefits of using three asymmetries, A^\q 2 ), A ( ^ a \q 2 ) and 
A^ e \q 2 ), built up from the amplitudes A±(q 2 ) and A\\(q 2 ) and interesting for studying the 
presence of physics BSM in the B — > K*£ + £~ decay. In fact this is the maximal number of 
independent asymmetries invariant under the transformations discussed in ref. [T7j. Involv- 
ing only Aj_(q 2 ) and A^(q 2 ), these asymmetries are less prone to hadronic uncertainties. 
We show that a study of their low (^-dependence (q 2 < rn 2 ^) can help discerning among 

various new physics scenarios, 
at low g 2, s and the point at which the three asymmetries can change the sign, one can tell 
which operators receive contributions from physics BSM. 

The asymmetry A^ m \q 2 ) has a pleasant property that it is exactly zero in the SM, and 
in the presence of a new physics phase its value can become significantly different from 
zero, and therefore its experimental study would be highly welcome. 

Furthermore, we introduced Aj, e \q 2 ), which is a simpler quantity than the commonly 
used A^\q 2 ), and it is related to it via 

24V) 

Since our quantities involve only A± ,\\(q 2 ), we could extend the expression for 4?V), 
derived in ref. [17], and write A^\q 2 ) in terms of functions h{q 2 ) entering the angular 
distributions of B — > K*£ + £~ to the massive lepton well. If the new physics does 

not alter the sign of Re Cj or Re Cg, we note that our asymmetry Aj, e \q 2 ) reaches the 

11 ^4^? \q 2 ) and other transverse asymmetries have been discussed in the literature [31 El [TTl [22] . Here we 
focus onto those that we consider phenomenologically more interesting because their smaller sensitivity to 
hadronic uncertainties. 



In particular from the shapes of A^\q 2 ) and A^^^q 2 ) 



+ 



4V 



+ 



A^ m \q 2 



+ 



Ap\q 2 
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maximum at low q 2 and is equal to one in the SM, while it gets suppressed in the presence 
of right handed currents and/or the new physics phase. 

Measuring all three asymmetries discussed in this paper should be within reach at 
LHCb and the Super-B factories. 
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Appendix 



In the numerical analysis we used m B = 5.279 GeV, and m K * = 0.892 GeV. For the SM 
Wilson coefficients we take [3]: 



Ci = -0.257, C 2 = 1.009, C 3 = -0.005, C 4 = -0.078, C 5 = 0.000 
C 6 = 0.001, Cf = -0.304, Cf - Y{q 2 ) = 4.211, C$j = -4.103, 



(56) 



computed to next-to-next-to leading logarithmic accuracy in MS(NDR) renormalization 
scheme at the scale /i = 4.8 GeV [2J. For the b and the charm quark mass we take [T9"] : 



m c MS (m c ) 



1.28(8) GeV. 



m 



MS, 



m 6 )=4.19(^) GeV. 



(57) 



Concerning the form factors in the range of low q s, the ratio i? and the slope z, defined 
in eqs. (23) and (25) respectively, have been uniformly (not Gaussianly) varied within the 
ranges [I31IH1 EES]: 

R E (0.17, 0.23) GeV" 1 , z E -(0.027, 0.033) GeV -2 . (58) 
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